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Abstract 

Drift-reduced plasma fluid models are commonly used in plasma physics for analytics and simulations; 
however, the validity of such models must be verified for the regions of parameter space in which tokamak 
plasmas exist. By looking at the linear behaviour of drift-reduced and full-velocity models one can 
determine that the physics lost through the simplification that the drift-reduction provides is important 
in the core region of the tokamak. It is more acceptable for the edge-region but one must determine 
specifically for a given simulation if such a model is appropriate. 


1 Introduction 

Fluid models have been used to describe plasma behaviour in a magnetic field since Braginskii derived the 
plasma fluid equations and calculated the transport coefficients in 1965 [T|. Many other fluid systems have 
been derived based on these original equations using various simplifications to describe any specific physics 
of interest. One such simplification that is often used is called the drift-reduction, in which the momentum 
equation is reduced by taking its curl resulting in an equation for the evolution of vorticity, W = V x v. 
For this new system to be closed, an assumption is made that the perpendicular velocities are dominated 
by the E x B drift, which relates the parallel vorticity to the potential: W\\ = (cgs Gaussian units 

are used for the duration of the paper). 

When this technique is used, time-scales below the ion cyclotron time are averaged over so fast waves 
are removed from the system. Also, the effect of the pressure gradient on the velocity evolution is not 
carried into the vorticity equation as the curl of a gradient is zero. This is especially significant because the 
largest pressure gradients will be perpendicular to the magnetic held and it is the perpendicular velocity 
equations which are replaced by an equation for the vorticity. 

These dropped physical effects will only be important in particular regions of parameter space, so this 
reduction is valid only when the missing physics is negligible. By using linearisation techniques, the be¬ 
haviour of these models can be compared to determine in which cases the drift-reduction is acceptable. 
The systems will be simplified to the incompressible limit (such that V ■ v = 0) to look at the most basic 
case that will still produce drift-waves. Any differences for this case, then, are fundamental and will carry 
on into more complex scenarios. Tokamaks are operated in well-defined yet broad parameter spaces, so 
the application of drift-reduced plasma fluid models for tokamak modelling can to be explicitly explored. 


2 Full velocity vs drift-reduced models 

A full-velocity model is one that evolves all three components of equation [lj the ion momentum equation 
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where m is the mass of the ions, n is the ion density, v is the ion velocity, W is the ion vorticity, J is the 
current density, B is the magnetic field, p is the total pressure, Xv is the velocity diffusivity, and c is the 
speed of light. By taking the curl of equation [lj an equation for vorticity is obtained. It is convenient 
to take the parallel component of the vorticity equation, as shown in equation [2j because it includes the 
behaviour of the perpendicular velocities: 
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where x and z are the perpendicular directions and y is parallel to the magnetic field line. This geometry 
will be explained in more detail in the next section. The normalised drift reduced equations, as derived by 
Hazeltine, et al. [2|, can then be written in the incompressible limit as follows: 
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where p is the pressure, is the electric potential, An is the parallel vector potential, [/, g] = — 

are the standard advection brackets, Jy = Vj_Am, W\\ = and rj is the parallel resistivity. The equation 

for U|| does not couple to these in the incompressible limit so is omitted, but it is important to note that 
the parallel velocity will evolve to maintain V ■ v = 0. 


Using the same normalisations, the full-velocity model in the incompressible limit is given by 
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with V • J = 0, J = V x B, and B = V x A. To evolve the vector potential, the generalised Ohm’s law 
given by Lifshitz j3] is used with temperature gradients neglected due to the isothermal assumptions. The 
equation for A x is excluded above because it simply evolves to maintain force balance without coupling 
to the remaining equations. For both systems parallel derivatives are taken along the perturbed field by 


defining Vy / = dy / - [A \\, /]. 


It is important when exploring the effects of the drift-reduction to ensure the two systems (full-velocity 
and drift-reduced) are identical in all other aspects. To do this, the full-velocity system was drift-reduced 
and in the linear limit reproduces exactly the dispersion relation of the Hazeltine model, which can be seen 
in the next section in equation [5j In this way, the effects of the gyro-viscous cancellation which is used in 
both models are not observed in our comparison. 
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Figure 1: The geometry for the linearisation is quasi-3D with equilibrium pressure gradient, current density, and 
magnetic field that satisfy force balance. Perturbations are in y and z such that the total perturbation is at an angle 
to the magnetic field, B 0 - 


3 Linearisation 

For all linearisations, we define a quasi-3D, orthogonal coordinate system ( x-y-z ) such that the equilibrium 
magnetic field Bq is in the y-direction, the equilibrium current density Jo is in the negative ^-direction, 
and the background pressure gradient is in the x-direction; however, perturbations are only in y and z 
with no extent in x consistent with a local approach, as detailed in figure [l] 

The background pressure gradient is present to drive the drift-wave instability, and the background cur¬ 
rent density and magnetic held are provided to satisfy force balance. All perturbations are of the form 
/ _ e i(k y y+k z z-fit) w ] iere q i s the complex frequency defined as ft = u + ij. Both systems were un¬ 
normalised prior to linearisation so that physical parameters, such as Alven speed, could be more easily 
substituted into the resulting dispersion relations. 


3.1 Drift-reduced dispersion relation 

The drift-reduced system in equation [3j once linearised, results in the following dispersion relation: 
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where w* is the drift-wave frequency, uo c i is the ion cyclotron frequency, ui p i is the ion plasma frequency, 
and va is the Alfven speed. These are defined as 
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Parallel Alfven waves as well as resistive drift-waves can be seen in the terms of the dispersion relation. In 
the case of zero resistivity the waves are stable and simply propagate. For tj > 0 the most unstable growth 
rate and corresponding frequency can be extracted using typical values for magnetic held, pressure gradi¬ 
ent, and background density within a tokamak (B =1T, 'S/p =10 6 Pa/m, no =10 18 m“ 3 , and k y = k z = 1). 
These values are also chosen to satisfy w* < u> c i such that the ion cyclotron frequency is the highest fre- 
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Figure 2: Growth rate (dashed) and frequency (solid) of the drift-wave instability as a function of conductivity for 
the drift-reduced system. Frequency, growth rate, and conductivity are all normalised to the ion cyclotron frequency. 


quency in the system (note that the pressure gradient here is typical for the pedestal and will be lower in 
other areas of the tokamak, reinforcing this ordering). The resulting frequency and growth rate are plotted 
in figure [ 2 ] as a function of conductivity, a = \jr\ (ie. the inverse of resistivity). 


3.2 Full-velocity dispersion relation 

When the full-velocity system is linearised, the resulting dispersion relation (equation [6]) contains the extra, 
fast-physics that was lost in the drift-reduction. 
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The perpendicular resistivity has been approximated to be double the parallel resistivity, as given in Wesson 
[5]. Notice this expression has a higher order in 17 compared with equation [ 5 ] due to the additional equation 
for the perpendicular velocities, which results in an extra mode in the growth rate and frequency. The two 
dispersion relations, equations [ 5 ] and [6j are identical when 17 4 term, the last three terms in 17 3 , and the 
parallel wave number in the second term of the Q 2 are neglected, indicating that these terms contain the 
physics lost in the drift-reduction. This includes various propagating parallel and perpendicular resistive 
modes and the ion cyclotron wave. 


For the same values of magnetic field, background density, and pressure gradient the growth rate and 
frequency of this expression are quite similar to that of the drift-reduced model, so it was necessary to look 
at the solutions over a large parameter space in conductivity and electron beta, defined by 

n _ Pgas _ &irnT e 

Pe ~ Pmag ~ ' 

In figure [3] the magnetic field is set to constant B = IT and the density is adjusted to vary beta. This 
is useful to do because the terms in equation [6] are not functions of only f3 e - they depend on various 
combinations of density and magnetic field. In essence, the parameter space is 3D, however this is not 
easily visualised so magnetic field has been held constant for illustrative purposes. 
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Figure 3: Full velocity and drift-reduced growth rates and frequencies as a function of conductivity and f3 e at 
Bq = IT. Growth rates and frequencies are normalised to the drift-wave frequency, w*. 


4 Tokamak relevance 

The parameter space in which tokamaks operate is specific to the region within the tokamak (core vs edge) 
and the particular tokamak in question. For a large tokamak of size similar to JET, the Joint European 
Torus at the Culharn Science Centre, the core operates around /3 e = 0.03 and a = 10 15 s _1 , while in the 
edge f3 e = 0.005 and a = 10 12 s _1 . 


Figure [4] depicts the percent difference in growth rates between the full-velocity and drift-reduced sys¬ 
tems as given by 
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It is clear that at low conductivity the drift-reduction breaks down for all values of j3 e . This is due to the 
terms exclusively in equation [g] that are functions of r/ and rj 2 becoming very large at low conductivity, 
a = rj~ l . 


At low /3 e , which corresponds to low density, we see a fairly universal disagreement between the mod¬ 
els. Since u* oc n -1 and va oc n -1 / 2 these plasma parameters become larger at low density. The plasma 
frequency u) p i oc n 1 / 2 , so it becomes small at low density, while the ion cyclotron frequency is not a function 
of density. All of the additional terms in equation [6] vanish at low density except for the second term of the 
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D 3 term, which is proportional to n 1 . At low density and high conductivity this term dominates, 

^ci 

but as conductivity is lowered, the tj 2 term takes over, thus the small area of agreement even at low f3 e . 



Figure 4: The percent difference from equation [ t] as a function of conductivity and electron beta. Usual operational 
regimes for the core and edge in JET are marked. 


The indication is that drift-reduced models are able to accurately reproduce edge behaviour, where errors 
can be as low as 0%. There are regions even in the edge, however, where a drift-reduced model may not be 
appropriate and errors can reach as high as 100%. In the core, there is a fairly consistent error of around 
100 % from the full-velocity solution, so full-velocity fluid or gyro-kinetic models should be used in this 
region. 


5 Conclusion 

Drift-reduced models provide simplified dispersion relations for more succinct analytics and the exclusion 
of fast waves allows for larger timesteps leading to faster simulations, so these models are an important 
subset of the full fluid description. The validity of these models has been tested for a simple quasi-3D slab 
resulting in drift-wave growth rates and frequencies that only agree with the full-velocity fluid description 
in specific regions of parameter space. Though the worst agreement lies outside of the operational regime 
of tokamaks, there is still questionable agreement for core simulations. It is important when choosing a 
model to use for tokamak plasma simulations to identify the parameter space in which the simulation will 
be operating as to identify whether a drift-reduced model is appropriate or if a more accurate, full-velocity 
model should be used instead. 
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